Investigation of mechanical losses of thin silicon 
flexures at low temperatures 

R. Nawrodt 1 ' 2 , C. Schwarz 2 , S. Kroker 2 ' 3 , I. W. Martin 1 , F. 

Briickner 3 , L. Cunningham 1 , V. Grofie 2 , A. Grib 4 , D. Heinert 2 , 

J. Hough 1 , T. Kasebier 3 , E. B. Kley 3 , R. Neubert 2 , S. Reid 1 , S. 

Rowan 1 , P. Seidel 2 , M. Thiirk 2 , A. Tiinnermann 3,5 

1 Institute for Gravitational Research, University of Glasgow, G12 8QQ, United 
Kingdom 

2 Institut fur Festkorpcrphysik, Friedrich-Schiller-Universitat Jena, Hclmholtzweg 5, 
D-07743 Jena, Germany 

3 Institut fur Angewandte Pliysik, Friedrich-Schiller-Universitat Jena, 
Albert-Einstein-Strafie 15, D-07745 Jena, Germany 

4 Kharkiv V. N. Karazin National University, Physics Department, 61077 Kharkiv, 
Ukraine 

5 Fraunhofer-Institut fur Angewandte Optik und Feinmechanik IOF, 
Albert-Einstein-Strafic 7, D-07745 Jena, Germany 

E-mail: r . nawrodt@physi.es .gla. ac .uk 
Abstract. 

The investigation of the mechanical loss of different silicon flexures in a temperature 
region from 5 to 300 K is presented. The flexures have been prepared by different 
fabrication techniques. A lowest mechanical loss of 3 x 10 -8 was observed for a 130 fim 
thick flexure at around 10 K. While the mechanical loss follows the thermoelastic 
predictions down to 50 K a difference can be observed at lower temperatures for 
different surface treatments. This surface loss will be limiting for all applications using 
silicon based oscillators at low temperatures. The extraction of a surface loss parameter 
using different results from our measurements and other references is presented. We 
focused on structures that are relevant for gravitational wave detectors. The surface 
loss parameter a s = 0.5 pm was obtained. This reveals that the surface loss of silicon 
is significantly lower than the surface loss of fused silica. 



PACS numbers: 00.00, 20.00, 42.10 



1. Introduction 

Gravitational wave (GW) detectors currently under operation are amongst the most 
sensitive instruments ever built sensitive for displacements. They are measuring tiny 
fluctuations of the space-time based on an interferometric principle. The optical 
components of a Michelson-like interferometer are suspended as pendula in order to 
reduce seismic coupling. Current detectors like LIGO p], Virgo [2] or TAMA [3] use 
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metal wire to suspend the test masses. GEO600 [I] pioneered the use of a monolithic 
suspension system based on fused silica fibres [5] Ej and hydroxide catalysis bonding 
[3 IB]. Second generation detectors like Advanced LIGO [9] or Advanced Virgo [10] 
will adopt this technique among others and aim for an increase of sensitivity of about 
an order of magnitude compared to initial detectors. Currently, a design study for a 
3 rd generation detector is underway in Europe [H] . In order to further increase the 
sensitivity beyond the 2 nd generation detectors a dramatic change in size, topology and 
materials for the detector are needed (see e.g. [12] ). Silicon is a promising candidate 
material for such a GW observatory due to its excellent thermal and mechanical 
properties [131 E]. Silicon would allow the fabrication of large optical substrates as 
well as suspension elements in a quasi-monolithic suspension [151 IE]- Furthermore, first 
attempts have been made to demonstrate a full reflective coating based on monolithic 
waveguides [T7] . 

The mechanical loss of a material is a critical parameter to estimate the thermal 
noise performance of an optical component or a suspension (see e.g. [18J ) . However, 
crystalline materials have in general a larger coefficient of thermal expansion compared 
to amorphous materials like fused silica. This results in a large thermoelastic noise at 
room temperature [19]. The use of crystalline materials like silicon makes it necessary 
to operate the detector at low temperatures to reduce the contribution of thermoelastic 
noise. 

It has been shown in the past that the mechanical loss of small structures can be 
deconvoluted into a bulk material intrinsic loss and a surface loss [201 EEJ [22]. This 
surface loss can be a significant and limiting source of mechanical loss - especially 
in small scale structures like suspension elements or the micro-structured surfaces of 
optical elements like monolithic waveguides. Our investigation focuses on structures 
which sizes that are relevant for the construction of suspensions and optical micro and 
nano structures. Furthermore, we exclusively focused on surface preparation techniques 
that are applicable for a long time operation in vacuum. 

2. Sample preparation and experimental technique 

The mechanical loss of silicon flexures was measured using a ring-down technique 
described elsewhere in detail [231 El]- The flexures were fabricated from high purity 
Si(100) wafers with their long axis parallel to the Si(110) direction. All flexures have a 
thicker block (525 /im) on one side and a thin flexure (75-130 /im) with lengths between 
35 and 50 mm. The thin flexures were obtained by two techniques: wet chemical etching 
and dry chemical etching. When the thin flexure is oscillating the thicker block at one 
end will only contain a reduced amount of vibrational energy allowing the flexures to 
be clamped at this side without affecting the mechanical loss significantly. Tabled] 
summarizes the parameters of the samples used in this paper. 

Due to the different fabrication techniques the surface topology of the samples is 
also different (see figure[l|). The wet chemical etching shows etch pits. The dry etched 
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Table 1. Summary of the parameters of the Si flexures. 



Properties 


Sample 1 


Sample 2 


Sample 3 


Sample 4 


length (mm) 


35 ±0.2 


50 ±0.2 


50 ±0.2 


50 ±0.2 


width (mm) 


4±0.2 


8 ±0.3 


8 ±0.3 


8 ±0.3 


thickness (fim) 










flexure 


75 ±2 


100 ±5 


130 ±10 


130 ±5 


block 


525 


525 


525 


525 


etching 


wet 


dry 


dry 


dry 


roughness 










front 


33 nm 


6.8 nm 


5.9 nm 


6.2 nm 


back 


0.9 nm 


1.9 nm 


1.6 nm 


337 nm 





(a) 



(b) 



Figure 1. Surface profile of a wet etched (a) and a dry etched (b) sample. Both 
images have been obtained with the same magnification and show an area of 600 fim 
x 450 jim. 



samples have a smooth surface with small (fim scale) droplets. The roughness of the dry 
etched samples is about a factor of 5 smaller than for the wet chemical etched. Sample 
4 was prepared from a single side polished wafer. The etching process was carried out 
from the polished side. The result was a silicon flexure with an identical geometry to 
sample 3 but one unpolished side. This sample was used to investigate the influence of 
the surface on the measured mechanical loss. 

Internal resonant modes of the samples were excited by means of an electrostatic 
drive plate and an applied high voltage (up to 1600 V). The excitation using an 
appropriate driving frequency allows the selection of different modes. The resonant 
frequencies /„, of the flexures can be calculated in the case of bending modes by means 
of ESI: 



fn 



Y 



1) 



2vrL 2 V 12p 

t is the thickness of the flexure, L its length, Y the Young's modulus and p the density 
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of the material. a n is a numerical factor and can be obtained from: 

cosa n x cosha n = —1 (2) 

The first four values are 1.875, 4.694, 7.855 and 10.996. 

The vibration of the silicon flexure was observed by a simple readout which senses 
the motion of a laser beam reflected from the surface. When a sufficiently large 
amplitude was reached the excitation was switched off and the subsequent ring down 
was recorded. The characteristic ring down time r (which corresponds to the amplitude 
decay to 1/e of the initial value) leads immediately to the mechanical loss 4>: 

^ = vr, (3) 

with the resonant frequency fo of the mode under investigation. The whole setup 
was placed into a specially built cryostat [26] providing the desired temperature and its 
stability over the measuring period. At low temperatures a typical temperature stability 
of better than 0.3 K was achieved within the ring down period. Typical ring down times 
r in the order of several 10 up to 1000 s have been observed for frequencies between 1 
and 86 kHz. All modes have been measured several times during 3 temperature cycles 
from 5 to 300 K. Each sample was re-clamped twice to reduce systematic errors arising 
from bad clamping. 

3. Mechanical loss in silicon flexures 

The measured loss of a silicon flexure will consist of different types of dissipation 
that causes energy loss. Being an integral measurement it is important to carefully 
study possible loss sources in order to be able to interpret the mechanical loss spectra 
throughout the full frequency and temperature region. The most likely types of 
mechanical dissipation occurring in the presented experiment are thermoelastic losses, 
phonon-phonon scattering losses, losses that are related to surface states and impurities 
or other losses that shall be summarized as excess losses. 

3.1. Thermoelastic loss 

Zener [27] was the first to investigate a loss mechanism due to dissipative heat flux in 
samples under bending oscillation. When a sample is bent parts are locally heated and 
cooled due to the inversion of the thermal expansion law. These local temperature 
differences cause heat fluxes that are accompanied by an increase of entropy - or 
dissipation. The temperature distribution T in an object can be obtained from the 
stress field cr^ by means of the generalized equation of heat conductivity [28], [29] : 
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and aiij are the tensors of thermal conductivity and thermal expansion. T is the 
average temperature of the sample. C p is the heat capacity. The mechanical loss induced 
by thermoelastic damping is then given by: 

4>te = 2£ J a ij d ij lm{f}dV (5) 

E is the total vibrational energy stored in the test body, T = Te tult and a = ae lujt with 
oj = 2nf the angular frequency of the mode under consideration. Eq. (jlj) and eq. (JSJ can 
be solved by means of a finite element analysis for arbitrarily shaped test bodies and 
arbitrary mode shape. 

Eq. §5§ can be derived analytically in the special case of pure bending modes in 
isotropic materials. For these modes the restoring force for the oscillation is exclusively 
based on bending forces. The simplified equation is given by 

a 2 YT our 



(6) 



with 



(7) 



P C p 1+U 2 T 2 

pC p t 2 

T = 

71 fi 

where t is the thickness of the sample. 
3.2. Phonon-phonon loss 

Another possible loss source in silicon is the phonon-phonon loss. In thermal equilibrium 
the phonon distribution is defined by the temperature T of the sample. Applying an 
external oscillation with a typical wavelength much longer than the phonon wavelength 
(like throughout this paper) results in a modulation of the lattice. This modulation 
shifts the phonon distribution. The process of reassembling all phonons to this new local 
equilibrium generates entropy and is thus a loss mechanism which is called Akhiezer loss. 
The mechanical loss associated with this mechanism is given by |30j : 



_ C P TY urrph , s 

V 2 l+U 2 T 2 h 

with the parameter 7 = 2.2, the speed of sound v and the mean phonon lifetime r p h 
which can be obtained from the thermal conductivity [30]. This assumption is not 
correct at very low temperatures where the thermal conductivity is limited by the 
sample dimensions [21] • However, this treatment will allow a comparison of the different 
contributions. 
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3.3. Surface loss 

Several authors have pointed out that a mechanical loss contribution of a micro cantilever 
is dominated by a thin surface layer (see e.g. j2ni ETJ [22} [32j |33]). This area is assumed 
to have different mechanical parameters compared to the bulk values. The origin of this 
change in parameters can be manifold, like local lattice distortions, adsorbed materials 
at the surface, open bonds, surface roughness, etc, and is not fully understood. The 
surface loss contribution can be written as (see e.g. [20]): 



= 4>buik (l + Vyjg) ( 9 ) 

with the mechanical loss of the bulk material (ftbuik and the volume-to-surface ratio V/S. 
The dissipation depth d s can be written as: 

d s = - — ^7 — [ (j){n)Y{n)dn. (10) 



4»bulkYbulk JO 

Here, Ybuik is the bulk's Young's modulus, <fr(n) and Y(n) are the distributions of the 
mechanical loss and the Young's modulus within the thin surface layer, h is the thickness 
of the surface layer. This description already assumes that the inhomogeneity of the 
surface layer is only dependent on the depth n from the surface. 

The numerical factor /i in eq. ([9]) takes the geometry of the sample and the 
modeshape of the resonance under investigation into account. It is given by: 

!s _V 5Js£W?r 

e is the strain amplitude and V and S are the volume and the surface area. For a flexure 
with rectangular cross-section oscillating in a transverse vibration this leads to 



(12) 

1 + t/w 

with the thickness t and width w of the flexure. Thin flexures lead to a constant value 
of n = 3. In the general case eq. ( ITT]) can be solved numerically using a FEA program 
to estimate the strain amplitudes e. 

It is often convenient to write eq. ([9]) in a different form: 

S 

= ^bulk + Ct S ^y (13) 

with the surface loss parameter a s = <fibuikd s - 

A similar approach for the surface loss was made by Yasumura et al. [21] . Assuming 
a thin flexure they obtained for the mechanical Q-factor of a surface loss dominated 
sample: 

Qsurface T^T - 77 Qs- iX^) 

bd K 
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Q s is the mechanical Q-factor of the surface layer (which is the reciprocal of the 
mechanical loss), t the thickness of the flexure, 5 the thickness of the lossy surface 
layer and Y s the Young's modulus of the surface layer. Assuming that the surface layer 
has a similar Young's modulus as the bulk material and rewriting eq. (I14p in terms of 
losses leads to: 

6 

4> surf ace = S(j) s — . (15) 

This coincides with the more general eq. fl9]) in the case of thin flexures (fi = 3, V/S ~ 
2/t) and the limit that the surface losses dominate the bulk losses. The surface loss 
parameter a s can then be identified with the product <p s 5. 

3.4.. Excess loss 

All losses so far are directly linked to the sample being investigated. Additional loss 
mechanisms occur during the interaction of the sample with the measuring setup. These 
additional losses are called excess loss. This term is a summary of different types of loss 
like residual gas damping (see e.g. [351 [36]) or electrostatic losses from the driving plate 
[351 137] . These losses can be reduced by carefully setting up the experiment and have 
been checked before the measurements [23j . 

3.5. Summary of possible internal loss mechanisms 

Figure[2] summarizes the internal mechanical loss contributions for a typical flexure 
under investigation in this paper (sample 3 from table[T]) for a frequency of 10 kHz. 
The material properties have been obtained from standard databases [38] . 

Thermoelastic loss dominates at temperatures above approximately 50 K. At 
around 125 K the effect of thermoelastic loss is significantly reduced due to the vanishing 




Figure 2. Summary of possible intrinsic loss sources of a silicon flexure at 10 kHz. A 
surface loss of 2 x 10~ 8 was assumed for the plot. 
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coefficient of thermal expansion a of silicon at this temperature. Here, the phonon- 
phonon loss sets the achievable limit. 

Depending on the exact value of the surface loss it will dominate at low 
temperatures. In the intermediate temperature region between 25 and 50 K the Akhiezer 
damping is the dominating loss source. Thus, in a setup with neglectable excess loss it 
should be possible to investigate the surface loss at temperatures below 20 K (see section 

KM- 

4. Experimental results 

4-1. Mechanical loss at low temperatures 

For samples 1 to 3 several modes have been measured. Starting with the fundamental 
bending mode and using eq. (P) it is possible to find the bending mode frequencies. It 
was possible to follow the sequence of bending modes up to about 70 kHz. Additional 
modes resulted in a total number of resonant frequencies per run of more than 30. The 
higher frequency modes have only been measured at lower temperatures. Close to room 
temperature their mechanical losses are high and the corresponding ring down times too 
short to be measured with sufficient accuracy. Figure[3] (a) shows the dependence of the 
measured resonant frequency on the parameter a n . Plotting f n against a 2 n should follow 
a straight line as predicted by eq. (JTJ. The fit can be used to determine the thickness of 
the samples. 




temperature (K) 



(a) (b) 

Figure 3. (a) - Plot of the determined resonant frequency at 300 K against for the 
first 10 modes of sample 1. The linear fit represents the behavior predicted by eq. (TT]) 
indicating that only bending modes have been selected, (b) - Temperature dependency 
of a typical resonant mode (3 rd bending mode of sample 2). The temperature 
dependency is caused by the temperature dependent elastic constants. 

The resonant frequency is temperature dependent. This behavior is determined by 
the temperature dependence of the Young's modulus (more exactly: by the temperature 
dependence of the elastic constants). At temperatures below 50 K the resonant frequency 
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Figure 4. Experimental results of the mechanical loss of sample 1 (wet etching, size). 



is approximately temperature independent. All frequencies throughout this paper are 
thus given at low temperatures unless otherwise noted. 

FigureH] summarizes the results obtained for the sample 1. The loss of all modes 
is thermoelastically limited at temperatures above 150 K. Below that temperature the 
modes show an excess loss compared to the predicted thermoelastic limit. The 19.98 kHz 
mode shows an onset of the thermoelastic dip at 125 K and reaches the lowest loss at 
the lowest temperature of about 5 x 10~ 8 . All other modes had a lowest loss well 
above 10~ 7 . A loss peak occured in the results of the 19.98 kHz mode. This peak 
was just observed in this specific mode. However, all other modes have shown higher 
losses which might have covered the loss peak. Alternatively, this peak could have its 
origin from a resonant coupling to internal modes of the clamping structure. Figure[5] 
compiles the results for sample 2. In contrast to sample 1 this flexure was fabricated 
by means of dry etching. The modes follow the thermoelastic limit down to 150 K. 
Below that temperature they show evidence of the vanishing thermoelastic damping at 
around 125 K. The 3 rd bending mode at 1958 Hz shows an excess loss at temperatures 
below 125 K. Here, a coupling between the sample's resonant frequency and the clamp is 
causing the additional damping. The lowest mechanical loss of 4.2 x 10~ 8 was observed 
for the 18.9 kHz mode again at the lowest temperature. The mode with the highest 
observed frequency showed some deviations from the predicted thermoelastic loss below 
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Figure 5. Experimental results of the mechanical loss of sample 2 (dry etching, size). 



125 K. The reason for this behavior is the deviation from the pure bending which is the 
prerequisite for the validity of eq. fl6]). Higher modes show complicated mode shapes 
including more and more torsional contributions. These contributions do not cause 
a time varying change in the stress tensor. Pure torsion does not result in a volume 
change and thus does not create a dissipative heat flux (see eq. (jlj)). Thus, the amount of 
thermoleastic damping is smaller than predicted by the simplified equation. Comparing 
the results to those of sample 3 (figure[6]) it is obvious that the excess loss seems to be 
only present at low frequency modes. This was a general observation throughout our 
investigations. All low frequency modes showed a slightly higher unexpected mechanical 
loss. Results for low frequency modes have been published previously on wet etched 
samples [23]. The results for sample 3 reveal an even lower mechanical loss for this 
sample reaching 3 x 10~ 8 at around 10 K. 

The results suggest that thin samples have slightly higher losses than thicker. This 
is a first indication of surface dependent losses which predict an increasing loss for 
thinner samples due to the increased surface-to-volume ratio. However, it needs to be 
first clarified that the low temperature loss is caused by the surface quality. Therefore, 
two identical samples with different surface preparation have been measured (sample 3 
and 4). While sample 3 had a dry etched and a mechanically polished surface sample 
4 consisted of a dry etched and an unpolished (lapped) surface. Sample 3 had two 
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surfaces with a small surface roughness (see tabled]). In contrast, sample 4 had one 
surface with a 50 times larger roughness. The results are summarized in figure[7](a). 
A log- log-plot was chosen to emphasize the low temperature part of the measurement. 
Both samples show a similar behavior for temperatures above 50 K. They both follow 
the thermoelastic limitation. Below 50 K sample 4 has a lowest mechanical loss above 
10 -7 whereas sample 3 showed a minimum loss of 3 x 10 -8 . Both measurements have 
been obtained at similar conditions: 10 K and the 18.4kHz mode (6 th bending mode). 
This measurement is evidence that the surface of the sample, rather than the setup, is 
limiting the mechanical loss at low temperatures. 

The wet etched sample had a different geometry and thus it is not possible to 
compare the results at the same frequency and modeshape. However, it is possible to 
compare the results for the 19.98 kHz mode of the wet etched sample with the dry etched 
samples although it is not the same modeshape. Due to the observed loss peak with 
unknown origin (see figureH^c)) only the lowest measured losses at temperatures below 
14 K were plotted in figured^) for comparison. The obtained loss lies in between sample 
3 and 4. This correlates with the roughness of the wet etched sample of 33 nm which is 
in between the values for the other samples. The results suggest a correlation between 
the roughness of the sample and the mechanical loss obtained at low temperatures. 

4-2. Analysis of the surface loss 

In order to extract a value for the surface loss of silicon at low temperatures the 
surface loss model by Gretarsson and Harry (see eq. f fT3|) ) was applied to available data. 
Eq. ( 1T31 suggests a linear correlation between the mechanical loss and fiS/V. Figure[7(b) 
summarizes different mechanical losses obtained in several experiments. Only long 
time experiments at low temperatures have been included in this analysis due to their 
applicability to a use of silicon in GW detectors or other long time experiments. It is well 
known that a lower surface loss can be obtained for special treatments of the surface like 
heating up to 1000 degrees for several seconds [391 EES [33]. However, this lower surface 




Figure 6. Experimental results of the mechanical loss of sample 3 (dry etching, size). 
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Table 2. Summary of the parameters of the oscillators from figure[7][b) used to obtain 
the surface loss parameter from literature values. 



point 


T (K) 








geometry 


reference 


1 


4.8 


1.5 x 10" 


-4 


2.98 


220 [im x 5 fxm x 0.06 fim 


SB 


2 


4.2 


4.0 x 10" 


_ g 


2.99 


300 [im x 10 [im x 0.07 urn 


EP 


3 


4.2 


1.2 x 10" 


-5 


2.86 


260 [im x 3.9 \im x 0.29 urn 


M 


4 


6 


5.0 x 10- 


-6 


2.94 


470 /im x 45 /im x 1.5 /im 


m 


5 


7 


5.5 x 10" 


-8 


2.96 


sample 1 




6 


10 


5.0 x 10- 


-8 


2.96 


sample 2 




7 


10 


3.0 x 10" 


-8 


2.97 


sample 3 




A 


7 


1.2 x 10" 


-8 


2.98 


0100 mm x 0.5 mm 


m 


B 


5 


2.2 x 10" 


-9 


2.36 


076 mm x 12 mm 


SSj 


C 


2 


5.0 x 10- 


-10 


0.77 


0106 mm x 229 mm 


m 



loss disappears after some seconds or minutes due to a possible surface contamination 
that covers the surface after that time. Due to the expected operational time of hours, 
days or even months and years for long time experiments like the GW detection these 
values have not been included into our consideration. Additionally, only values from 
the literature have been chosen where it was possible to extract the geometry and the 
modeshape to calculate the geometry factor \i. Only experiments have been selected 
on silicon flexures at low temperatures (T < 10 K) to be able to neglect thermoelastic 
damping. The numbers at the plot indicate the original reference for the values. All 
key parameters of the different experiments are compiled in table[2j 

Plotting the mechanical loss of a silicon flexure based oscillator against its volume- 




Figure 7. (a) Investigation of the influence of the surface roughness on the mechanical 
loss. The polished sample corresponds to sample 3 and the unpolished to sample 4 
from table [TJ Both plots are obtained from the 18.4kHz mode (6 th bending mode), 
(b) Summary of the lowest loss obtained from different sized silicon oscillators. The 
numbers and letters are explained in tableland the text. The plotted line corresponds 
to a surface loss parameter a s of 0.5 pm. 
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to-surface ratio V/ S results in a linear dependency if the mechanical loss is determined 
by surface loss. The mode shape dependent factor fi was obtained from eq. ( 1 1 2 p for all 
Si flexures (1-7). The geometry factor for the bulk samples (A-C) has been obtained 
from a numerical derivation of eq. ( TTTj) using the FEA software COMSOL. \i is close to 
3 for all cantilevers and the thinnest bulk sample. The thick samples show a deviation 
which is largest for sample C. Here, the thickness is much larger than the diameter. 
This causes only a small fraction of the test body vibration to probe the surface loss. 
The effect of the surface loss to this specific mode is thus reduced. 

Eq. (flUl) has been used to model the data obtained for the thin flexures (1-7). It 
was assumed that the intrinsic bulk loss is much smaller than the surface loss. The only 
free parameter in eq. ffl3|) is the surface loss parameter a s under this assumption. It 
was adjusted in a way that the resulting line corresponds well to the lowest losses of the 
silicon flexures. This allows a rough estimate of the surface loss parameter of 0.5 pm 
with an error of about 25%. 

The mechanical loss values obtained from the flexures follow the predictions of the 
surface loss model. The bulk samples show a deviation towards higher losses. Here, the 
neglected intrinsic bulk loss starts to become significant. It is possible to estimate an 
averaged surface loss 4> sur f assuming a homogeneous surface layer of thickness t s and a 
Young's modulus equal to the bulk value. Eq. (TTOj) gives then: 

d s = t^L. (16) 

(pbulk 

This leads to a value of 

surf ace C^s/^s (17) 

for the surface loss. The thickness of the lossy layer t s is unknown and can only be 
approximated. Assuming similar values compared to silica (t s « 1 /im, [20]) gives a 
homogeneous surface loss value of 5 x 10 -7 which is a factor of 20 smaller than for fused 
silica [20] . However, this estimate is based on the weak assumption of a similar thickness 
of the lossy surface layer in silicon and silica. 

5. Conclusions 

We presented measurements of the mechanical loss of silicon flexures in a temperature 
range from 5 to 300 K and frequencies from 1 to 86 kHz. At temperatures above 100 K the 
experimental data follows the thermoelastic predictions. The observed loss is higher than 
the thermoelastic predictions for temperatures below 50 K. The level of the mechanical 
loss at low temperatures is strongly dependent on the surface quality. Rough surfaces 
produce a higher mechanical loss than etched smooth surfaces. The lowest mechanical 
loss of 3 x 10 -8 was obtained for a 130 /im thick cantilever at a frequency of 18.4 kHz 
and a temperature of 10 K. Using the lowest mechanical loss obtained from different 
silicon based oscillators it was possible to extract a surface loss parameter according 
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to the surface loss model from Gretarsson and Harry of a s = 0.5pm. This value is at 
least an order of magnitude smaller than for fused silica. Together with its remarkable 
mechanical and thermal properties the small surface loss makes silicon a very promising 
material for low mechanical loss applications at cryogenics like suspension elements in 
future gravitational wave detectors. However, the origin of the surface loss is still not 
fully understood and will be within the focus of further investigations. 
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